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ABSTRACT 

We present high resolution, long-slit spectra of the jet-cloud interaction in the powerful 
southern radio galaxy PKS2250-41. We have resolved the emission lines into two main 
kinematic components: a broad component (FWHM>900 km s~^) and a narrow com- 
ponent (FWHM<150 km s^^). While the broad component is characterized by a low 
ionization level (with particularly weak HeIIA4686 emission) and is spatially associ- 
ated with the radio lobe, the narrow component is characterized by a higher ionization 
level and extends well beyond the radio lobe. Crucially, we measure a higher electron 
temperature for the broad component (T '^30,000 K) than for the narrow component 
(r ^15,000 K). The general line ratios and physical conditions of the two components 
are consistent with a model in which the broad component represents gas cooling be- 
hind the shock front driven by the radio jets, while the narrow component represents 
the AGN- or shock-photoionized precursor gas. However, uncertainties remain about 
the gas acceleration mechanism behind the shock front: unless the radio components 
are expanding unusually fast in this source, it is likely that entrainment of the warm 
clouds in the hot post-shock wind or radio plasma is required in addition to the initial 
acceleration across the shock front, in order to explain the large line widths of the 
broad component. 

The similarities between the kinematic properties of PKS2250-41 and some high 
redshift radio galaxies suggest that the ambient and the shocked gas have also been 
resolved in the more distant objects. Given the evidence that the emission line pro- 
cesses are affected by the interactions between the radio and the optical structures, 
care must be taken when interpreting the UV spectra of high redshift radio galaxies. 
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1 INTRODUCTION 

The study of Seyfert galaxies and low redshift radio galax- 
ies has revealed extended emission line properties which are 
well explained in terms of an illumination model {e.g. Fos- 
bury 1989) in which the ambient gas is illuminated and 
photoionized by active galactic nuclei (AGN) in the cores 
of the galaxies. This scenario has the advantage of being in 
agreement with the unification models {e.g. Barthel 1989) 
in which it is proposed that many active galaxies are in- 
trinsically identical but appear different due to orientation 
effects. 

However, AGN illumination is not the whole story. HST 
images of Seyfert galaxies (Capetti et al. 1996) have shown 
a close association between the narrow-line region emis- 



sion line structures and the radio structures. The emission 
line profiles also reveal disturbed kinematics and high velo- 
city motions {e.g. Axon et al. 1998). These properties show 
that strong interactions are taking place between the radio- 
emitting plasma and the emission-line gas. 

Some nearby radio galaxies and the majority of high 
redshift radio galaxies (HzRG) {z >0.7) present clear evi- 
dence for similar interactions. Many HzRG show highly colli- 
mated UV continuum and emission line structures which are 
closely aligned with the radio axis (Chambers et al. 1987; 
McCarthy et al. 1987; Best, Longair & Rottgering 1996). 
The emission line spectra show a highly disturbed kine- 
matics which cannot be explained in terms of purely gravita- 
tional motions {e.g. McCarthy et al. 1996, Rottgering et al. 
1997) but suggest that the gas is interacting with the radio- 
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emitting components. Recently, work by the Stromlo group 
has provided a theoretical framework for understanding the 
emission line spectra produced by warm emission line clouds 
involved in jet-cloud interactionsT] {e.g. Sutherland, Bicknell 
& Dopita 1993; Dopita & Sutherland 1996). 

Jet-cloud interactions are therefore a common phe- 
nomenon in active galaxies and can have a strong influence 
on the properties we observe. The way the interactions take 
place is, however, poorly understood and many questions re- 
main open. In particular, the gas acceleration mechanism in 
jet-cloud interactions remains uncertain, and the extent to 
which the jet-induced shocks ionize the emission line clouds 
is not yet clear. 

We are carrying out a project whose main goal is the 
study of the interaction between the radio and optical struc- 
tures in radio galaxies. The very large distances of HzRG 
make this study difhcult due to the small angular sizes, the 
faintness of the objects, and the fact that optical observa- 
tions sample the UV rest frame, where the emission line 
physics are highly uncertain. Therefore we are concentrating 
on radio galaxies at intermediate redshifts which are much 
easier to study. 

PKS2250-41 (z=0.308) is one of the few relatively 
nearby radio galaxies with clear evidence for jet-cloud inter- 
actions. Due to its proximity and the strength of its interac- 
tion, PKS2250-41 is an excellent target to study the physics 
of this phenomenon in powerful radio galaxies. What we 
learn from PKS2250-41 will provide important clues about 
what is happening in more distant radio galaxies. The results 
will be useful to understand the role played by interactions 
between the radio- and optical-emitting structures in active 
galaxies in general. 



1.1 The jet-cloud interaction in PKS2250-41. 
Summary of previous results. 

Clark et al. (1997) and Tadhunter et al. (1994) showed that 
the interaction between the radio jet and the emission line 
gas in PKS2250-41 has the following consequences: 

• Morphological associations between the radio and op- 
tical structures: The dominant feature of the emission line 
images ([OIII]A5007) is an arc-shape structure located at 

~5.5 arcsec (~33 kpc)r| west of the nuclear continuum cen- 
troid. This arc circumscribes the Western radio lobe. 

• Complex kinematics with relatively broad low ioniza- 
tion lines. An anticorrelation is observed between line width 
and the ionization state in the Western arc. This effect may 
be explained, as in SN remnants (Greidanus & Strom 1992), 
as due to the fact that low ionization emission lines are pro- 
duced mainly in the shocked, compressed and accelerated 
gas, while the high emission lines are mainly emitted by the 
precursor gas. 

• High pressures measured in the western radio lobe, were 
the interaction is strongest. This suggests that the radio 

* Throught this paper we use "jet-cloud interaction" as a generic 
term to denote interactions between the radio-emitting compo- 
nents — including the radio jets, hotspots and lobes — and the 
warm emission line clouds. 

T Ho = 50 km s^^ Mpc^^, qo = 0.0 assumed throughout, giving 
an angular scale of 6.0 kpc arcsoc"^ 
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Figure 1. [OIII] narrow band image of PKS2250-41 (from Clark 
et al. 1997) with the slit positions ovorplotted. 



plasma and the line emitting gas are in approximate pressure 
equilibrium. 

The authors propose a scenario in which a direct colli- 
sion between the radio jet and a companion galaxy in the 
surrounding group is taking place. 

In this paper we present new long-slit spectra for 
PKS2250-41, which cover a larger spatial region than stu- 
died in the previous work, and contain information not only 
for the regions along the radio axis, but also perpendicular 
to it. The high spectral resolution of these data has allowed 
us to study in detail the properties and origin(s) of the diffe- 
rent kinematic components revealed in the strong emission 
lines. 



2 OBSERVATIONS 

The spectroscopic observations were carried out on the 
nights 1994 September 27/28 using the Royal Greenwich 
Observatory (RGO) spectrograph on the Anglo Australian 
Telescope. The detector was a Tek CCD with 226x1024 pi- 
xels of 27 /iin^, resulting in a spatial scale of 0.81 arcseconds 
per pixel. 

The data were taken with the slit approximately paral- 
lel (PA270) and perpendicular (PAO) to the radio axis us- 
ing both intermediate (3.0 - 3.3A FWHM) and high (1.4 - 
I.7A FWHM) spectral resolutions. We covered most of the 
projected area of the Western arc by obtaining spectra along 
three slit positions at PA270: along the radio axis, 2.5 arc 
sec South (2.5" S slit hereafter) and 2 arc sec North (2"N slit 
hereafter). The slit positions are presented in Fig. 1. A log 
of the spectroscopic observations is shown in Table 1. 

We did the same with the overlaps between the cross cut 
and the 2.5" S and 2"N slits. We detect (as expected) three 
kinematic components (two narrow-fone broad) at the radio 
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Table 1. Log of the spectroscopic observing run. 



Central A Grating texp 



AA 



Resol.(A) SeeingC) Slit(" 



Radio axis (PA270) 

4883 1200V 3600 4485-5300 1.7±0.1 0.6 1.5 

6250 600R 5400 5415-7040 3.0±0.3 1.0 1.5 

6560 1200V 5400 6165-6980 1.4±0.1 0.6 1.5 

8300 600R 6000 7520-9130 3.3±0.3 1.0 1.5 



6560 


1200V 


5400 


2.5 arc sec S 
6165-6980 


1.4±0.1 


0.6 


1.5 


6560 


1200V 


5400 


2 arc sec N 
6165-6980 


1.4±0.1 


0.6 


1.5 


6560 


1200V 


5400 


Cross cut (PAO) 
6165-6980 


1.4±0.1 


0.6 


1.5 



axis position. We have isolated the redder narrow component 
and plotted its properties in Fig. 5 (solid triangles). The 
good agreement with the cross cut values for the narrow 
component confirms that the blue component ts not detected 
N and S of the radio axis along PAO. 

The data reduction was carried out using standard 
methods provided in IRAF (a detailed description can be 
found in Villar-Martm et al. 1998). 



3 ANALYSIS AND RESULTS 

3.1 The fitting procedure 

Both IRAF and STARLINK (DIPSO) routines were used to 
measure the emission line fluxes, line widths and line cen- 
ters. For the blends, decomposition procedures were used in 
STARLINK (DIPSO), fitting several Gaussians at the ex- 
pected positions of the components. As the number of pos- 
sible mathematical solutions is large, we applied theoreti- 
cal constrains when necessary (like fixed ratios between line 
fluxes or separation in wavelength). In some cases, the in- 
formation obtained from the strong [OIII]AA5007,4959 lines 
was used to constrain the input parameters for the flts to 
other lines. We assumed, for instance, that all lines have the 
same number of kinematic components as those discovered 
in the [OIII] lines. Sometimes we also constrained the kine- 
matic components to have the same velocity width (taking 
into account instrumental broadening and wavelength de- 
pendence of velocity) as the [OIII] components. 

The [OIII]AA5007,4959 and H/3 lines were used to study 
the spatial variation of the properties of the different kine- 
matic components revealed by the fitting procedures. In or- 
der to do this, we extracted a 1-D spectrum for each spatial 
pixel across the Western arc the region of greatest interest to 
us, which shows the clearest signs of jet-cloud interactions. 
Those pixels for which there was a risk of contamination by 
the emission from the near-nuclear regions of the host galaxy 
were excluded. As an example of the fitting procedure, we 
present in Fig. 2 the results of the fits to the [OIII] lines 
for spatial pixels across the Western arc (slit aligned along 
PA270). 

Apart from the [OIII]AA5007,4959 and H/3 lines, most 
of the lines are faint, and the spatial information is difficult 



to obtain on a pixel-by-pixel basis for the individual com- 
ponents. In order to increase the signal to noise ratio and 
to obtain information on the individual kinematic compo- 
nents, we spatially integrated the emission from the bright- 
est pixels across the Western arc (using 1.5x3.2 arcsecond 
aperture, centered 5.6 arcseconds to the west of the nucleus 
along PA270) and fitted the fainter lines in the integrated 
spectrum. The strongest lines in each spectrum were used 
to constrain the input parameters for the fits to the fainter 
lines in the same spectrum. 

The result of this procedure is the isolation of the diffe- 
rent kinematic components in each line, characterized by 
flux, line width and central wavelength. It is, therefore, 
possible to study the kinematics and line ratios of each kine- 
matic component separately. 



3.2 Results from the PA270 slits 

In this section we analyse the information provided 
by the [OIII]AA5007,4959 and H/3 lines in the high resolu- 
tion spectra obtained along PA270 (approximately parallel 
to the radio axis). Our goal is to isolate the different gaseous 
components contributing to the line emission and study the 
spatial variation of the kinematic and ionization properties. 
The results for the perpendicular slit position (PAO) are pre- 
sented separately in section §3.3. 

The fits to the lines reveal complex kinematics with the 
presence of at least three spatially extended components: 
two narrow components (FWHM~60-200 km s~^) separated 
by ~140 km s~^ and a broad component (^500- 900 km s~^). 
We will distinguish the two narrow components by labeling 
them as the blue narrow component (component at shorter 
wavelength) and the red narrow component (component at 
longer wavelengths). 

We plot the spatial variation of the kinematic, flux and 
ionization properties of the narrow components in Fig. 3. 
The properties of the broad components are presented in 
Fig. 4. In order to highlight the possible connection with 
the radio structures, we have indicated in these diagrams 
the position of the radio hot spot and the outer edge of the 
radio lobe. 

The narrow components (Fig. 3) 



© 1999 RAS, MNRAS 000, §-?? 



4 Villar-Martin et al. 



10.0"W 



9.2"W 



8.4"W 



7.6"W ^: 




6.8"W 




^^A^^^^^^'OTy^^^^^ i/U >^^^^^ 




^^..n^f^tir fflfU^^rt^^ 'III ^.yjj'i.'^^fl^ 



A+l 




I^H^KJ, 



6.0"W 



5.2"W 



4.4"W 



Figure 2. The fitting procedure: Example of the spectral decomposition procedure done for PKS2250-41. The central panel shows 
the 2-D [OIII]A5007 spectra along the radio axis. Dipersion in A varies in the vertical direction and space in horizontal direction. The 
plots show the 1-D spectra created from every column. The emission lines of every spectrum were fitted with the necessary kinematic 
components. The smooth solid lines indicate the results of the fits. The inner and outer pixels in the arc have been indicated with arrows. 
The position of each spectrum relative to the nucleus of the galaxy is also indicated. 
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Fig. 3 shows that aU the narrow components have 
similar line widths across the full extent of the nebulosity 
(FWHM -60-200 km s"^). The velocity shift with respect 
to the nuclear emission is rather constant for each compo- 
nent. The blue narrow component — detected along the ra- 
dio axis and also in the detached emission line region ~4.5 
arcseconds to the NW of the nucleus (see Figure 3) — is 
noticeably displaced in velocity with respect to the other 
components (velocity shift: 80-140 km s~^). 

As revealed by the [OIII]/H/3 ratio, the ionization level 
of the narrow components shows little variation, both in 
absolute value and in spatial variation, across the arc. Due 
to the faintness of H/9 it was necessary to add all the pixels 
along the 2.5" S slit. The integrated value for this slit fits in 
with the variation defined by the other narrow components. 

The fiux curves (Fig. 3 bottom right) of the redder na- 
rrow component along the radio axis and the narrow com- 
ponent along the extreme slit positions are almost identical. 
The absolute fiux values also remain rather constant when 
we move North and South from the radio axis. The redder 
narrow component extends far beyond the hot spot (6.5 arc 
sec) and ~4.8 arc sec beyond the edge of the radio lobe. 

In contrast, the fiux curve of the narrow blueshifted 
component detected along the radio axis is different in the 
sense that this component disappears abruptly at the edge of 
the radio lobe. 

The broad components. (Fig. 4) 

The broad component presents a strong variation in its 
kinematic, ionization and fiux properties along the radio 
axis. The fiux peaks at the position of the hot spot, where 
its emission dominates over the narrow components. It is 
also at this position that the line is broadest (FWHM~900 
kms"^). 

At other positions along the radio axis and the extreme 
slits, the FWHM and fiux are rather constant. The velocity 
curves are very similar along the radio axis and the 2"N slit, 
with an S-shape, although spatially shifted (~1.6 arc sec). 

The ionization level of the broad component is low com- 
pared to the narrow components ([OIII]/H/3 ~1.5-4). As ob- 
served for the narrow component, the ionization level in- 
creases towards the galaxy along the radio axis. 

It was necessary to bin together several pixels to in- 
crease the S/N ratio along the 2.5" S and 2"N slits. The 
error is still too large for the measurement on the 2"N slit, 
but the result obtained from the 2.5" S spectrum is consis- 
tent with a very low ionization level of the broad component 
([OIII]/H/3 ~2.8), similar to the value measured along the 
radio axis. 



3.3 Results from the cross-cut slit (PAO). 

In this section we analyse the information provided by 
the [OIII] AA5007,4959 and H/3 lines in the intermediate reso- 
lution spectra obtained with the slit perpendicular to the 
radio axis. 

Two spatially extended components are detected: a 
broad component (FWIIM~500-900 km s~^) and a narrow 
component (FWHM ~200 km s"^). The resuhs for both 
components are plotted in Fig. 5. 

The low spectral resolution of this spectrum prevents us 
from resolving the two narrow components detected along 



the radio axis (see §3.2). Therefore, the information near and 
on the radio axis position is contaminated by the presence 
of the blue narrow component along the radio axis. In order 
to eliminate this effect, we have extracted a 1-D spectrum 
from the high resolution spectrum obtained along the radio 
axis and corresponding to the spatial overlap with the cross 
cut slit (see Fig. 1). 

The narrow component (Fig. 5) 

The narrow component is detected across the full extent 
(~ 8.9 arc sec) of the emission line nebulosity in the N-S 
direction along PAO. The FWHM (100-150 km s"^), velocity 
shift (150-200 km s'^) and ionization level ([OIII]/H/3=6.5- 
8) are rather constant across the slit (N-S direction). This 
is clear once we have corrected for the contamination of the 
blue narrow component at the radio axis position. 

The slight increase on the FWHM at the position of the 
radio axis, the V shape of the velocity curve, and possibly 
also the drop in the ionization level are due to the contami- 
nation by the blue component. The flux curve shows also a 
drop across two pixels at the radio axis position. The drop 
must be more pronounced, since the blue narrow component 
contaminates the emission in these two pixels. This suggests 
that the narrow component presents an edge brightening in 
the direction perpendicular to the radio axis. 



The broad component (Fig. 5) 

The broad component is spatially extended by 4.9 arc- 
seconds N to S across the radio axis. 

The width varies substantially across the arc in the N-S 
direction and peaks at the position of the radio axis. The 
width drops in the outer parts of the arc and is consis- 
tent with the values measured along the 2.5" S and 2"N slits 
(FWHM~500 km s~^). The flux has a maximum near the 
radio axis position. 

The ionization level of the broad component is remar- 
kably low ([OIII]/H/3=1.2-3) compared to the narrow com- 
ponent, as expected from our previous results (see §3.2) 

In summary, disturbed kinematics is apparent across 
the whole of the W arc. Multiple components are detected 
which present different spatial extensions, kinematic and 
ionization properties. A summary of the results is presented 
in Fig. 6. 



3.4 Additional information along the radio axis. 

We have spectra along the radio axis covering a wider spec- 
tral range, including numerous optical lines which can pro- 
vide additional information on the origin and physical pro- 
perties (density, temperature) of the kinematic components 
and their physical properties. 

After all the 2-D frames were aligned, we extracted 1- 
D spectra from the brightest spatial pixels of the Western 
arc (1.5x3.2 arcsecond aperture, centered 5.6 arcseconds to 
the west of the nucleus along PA270). We then isolated the 
broad and narrow components in all the interesting emission 
lines, using the fitting method described above. 

Table 2 presents the fiuxes of the main optical lines rela- 
tive to H/3 for the two components revealed by the fits. No- 
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Figure 3. The narrow line emission along the slits parallel to the radio axis: Spatial variation of the kinematic, flux and ionization 
properties of the narrow components detected across the Western arc along the slit positions parallel to the radio axis. The spatial 
shift has been calculated with respect to the continuum centroid of the galaxy. Panels: Top left: FWHM. Top right: Velocity shift with 
respect to the nuclear emission. Bottom left: Ionization level. Bottom right: Fluxes. Symbols: Radio axis red component -> solid 

triangles. Radio axis blue component -> open triangles. 2.5" S slit -> x symbols. 2"N slit -> open squares. The dashed and the dotted 
vertical lines mark the positions of the hot spot and the edge of the radio lobe respectively at the slit located along the radio axis. Note 
the similarities of the narrow line emission along the 2.5" S and 2"N slits, and the redder narrow component along the radio axis. 



tice the clear differentiation between the broad and the na- 
rrow components: as the previous results suggest, strong low 
ionization lines dominate the spectrum of the broad compo- 
nent and high ionization lines dominate the spectrum of the 
narrow component. 



Table 2. Line fluxes relative to H/3 of the broad and narrow com- 
ponents measured along the radio axis in the Western arc. H/3 flux 
is given in units of W m~^. The measurements were made from 
a spectrum extracted using 1.5x3.2 arcsecond aperture, centered 
5.6 arcseconds to the west of the nucleus along PA270. 



Particularly notable are the relative faintness of 
HeIIA4686 and the low values of the temperature-sensitive 
[OIII]A5007/[OIII]A4363 line ratio measured in the broad 
component (see below). Although these characteristics were 
also noted by Clark et al. (1997) on the basis of single Gaus- 
sian fits to low resolution spectra, the isolation of the di- 
fferent kinematic components shows that these are properties 
of the broad component. Thus, it is likely that the ionization 
minimum deduced from low resolution spectra at the posi- 
tion of the radio lobe by Clark et al. (1997) is a consequence 
of the increased dominance of the low ionization broad com- 
ponent at that location. The new results also suggest that 
the anticorrelation between the line width and ionization 
stage deduced from single Gaussian fits to the low resolu- 
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Flux(H/3) 



(2.7±0.1)xl0" 



(7.1±0.2)xl0" 



[OII]A3727 

H7 

[OIII]A4363 

HeIIA4686 

[OIII]A5007 

[NI]A5199 

[OI]A6300 

Ha 

[NII]A6583 

[SII]AA6716+6732 



7.2±0.3 

0.45±0.06 

<0.1 

0.22±0.03 

5.4±0.2 

0.24±0.04 

0.8±0.1 

3.3±0.2 

1.9±0.1 

2.2±0.3 



4.7±0.07 

0.32±0.05 

0.07±0.01 

<0.08 

1.35±0.06 

0.22±0.03 

1.10±0.08 

3.1±0.1 

2.63±0.09 

3.4±0.2 



© 1999 RAS, MNRAS 000, |-?? 



PKS2250-41 and the role of jet-cloud interactions in powerful radio galaxies. 7 



s - 



: X ^^ "" 1 ^ 



1000 



800 



600 



200 







y. 






> 



e 



B 



ee 



2 4 6 

K<--ShiftC 



10 12 

-->w 



2 4 6 

E<--Shift(' 






o 



2 4 6 8 10 12 

E<--ShiftC')-->W 



2 4 6 

E<--Shift(' 



300 



200 



100 



-100 



-200 



10 12 

-->w 




10 12 

-->w 



Figure 4. The broad line emission along the slits parallel to the radio axis: Radio axis -> open circles. 2.5"S slit -> x symbols. 2"N slit 
-> open squares. Note the maximum of the FWHM and the flux of the broad component at the position of the hot spot. Note also the 
low ionization level of this component. 



Table 3. Results of the fit to the [OIII]A4363 line. For the ra- 
dio axis the measurements were made from a spectrum extracted 
using a 1.5x3.2 arcsecond aperture centered 5.6 arcseconds west 
of the nucleus along PA270, whereas for the cross-cut a 4.0 xl.5 
arcsecond extraction aperture centered on the radio axis was used. 

[OIII]A5007/A4363 Narrow Broad 



Radio axis 


>52 


17±5 


Cross cut 


50±25 


10±4 


Added 


50±20 


13±4 



tion spectra is a consequence of the low ionization state of 
the broad component. 



3.4-1 The electron temperature 

We have used the [OIII]A4363 and [OIII]A5007 lines to mea- 
sure the electron temperature of the two main kinematic 
components. H7 and [OIII]A4363 were fitted together. We 
have done this analysis for three different spectra: 1) spec- 
trum of the arc extracted along the radio axis from those 
pixels where the broad component was brightest (same aper- 



ture as in §3.4) 2) spectrum extracted along the cross cut 
slit position also from the pixels where the broad compo- 
nent was brightest (1.5x4.0 arcsecond extraction aperture 
centered on the radio axis); and 3) adding these two spectra 
after correcting for any possible shift in wavelength. 

The best fit was obtained by constraining the velocity 
widths (obtained from the [QUI] A A 5007,4959 lines in the 
same spectrum) of all the kinematic components involved. 
We show in Fig. 7 the results of the three fits. 

Although the [OIII]A4363 line is faint, we are confident 
that the results are reliable because: 1) the fits to all spec- 
tra produce consistent results; and 2) the H7/H/3 ratios are 
consistent within the errors with Case B recombination. 

We have used the measured [OIII]A5007/A4363 ratios 
to calculate the electronic temperature of the gas using the 
formula in Osterbrock 1989 valid for low densities. The den- 
sities in the arc (see Clark ei al. and next section) do not ex- 
ceed a few hundred and the low density formula is valid. We 
present in Fig. 8 the function [OIII]A5007/A4363=f(r), with 
the measured values of the line ratio for the two kinematic 
components (for the three fits) plotted. The broad compo- 
nent lies in a region of the diagram where [OIII]A5007/A4363 
is highly sensitive to the temperature. For this reason the 
uncertainty in the temperature is large, in spite of the 
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Figure 5. Spatial variation of the kinematic, flux and ionization properties of the arc along the crosscut slit position. Both the narrow 
(open triangles) and the broad (open circles) components are plotted. The spatial shift is now calculated with respect to the radio axis 
position (dashed vertical line). Note change of symbols: Broad component -> Open circles. The solid triangles (see text) are measurements 
from the PA270 spectra extracted from the regions where they overlap with the PAG spectrum (see Fig. 1). 



relatively small errors in the [OIIIJA5007/A4363 ratio. It 
is clear, however, that the broad component has a higher 
electron temperature than the narrow component: we ob- 



tain, Te = 30,000 



+8,000 



K, Te = 60,000 



-1-100,000 



K, and 



-4,000 '^^J ^e — uu, uwu_25ooo 

Te = 41,000_gQQQ K for fits to the radio axis, cross-cut 
and co-added spectra respectively, while for the narrow com- 



ponent we obtain Te = 15, 000_ 



Note that, unlike the 



situation for the nuclear regions of active galaxies, the equiv- 
alent width of the [OIIIJA4363 line is large in the regions of 
interest. Consequently, systematic effects related to the sub- 
traction of the underlying stellar continuum are unlikely to 
be a significant source of error. 

We have also calculated the temperature using the to- 
tal fiux of the lines in the co-added spectrum (with no 
constraint on the properties of the lines). The result — 
[OIII]A5007/A4363= 28±8 — is consistent with an interme- 
diate temperature between the narrow and the broad com- 
ponent (see Fig. 8): Te=21000l|[JE;[5K. This supports the con- 
sistency of our results. Assuming that the narrow component 
is cooler, Te=21000 K is a lower limit on the temperature of 
the broad component. 



3.4.2 The densities 

We have attempted to use the [SII]AA6716,6731 doublets to 
calculate the density (Osterbrock 1989) of the broad and the 
narrow components. We show in Fig. 9 the results of the fit. 

We obtained '^"'^f ^^ =1.47±0.191 for the narrow 

6731 

component (consistent with low density limit), and 
- — g^gj — = 1.13 ±0.08 for the broad component. These re- 
sults provide tentative evidence that the broad component 
has a higher density {rie = (5.7^^g) x 10* m~^, assuming 
Te=30,000) than the narrow component, although data with 
a higher signal-to-noise ratio and a better sky subtraction 
will be required to confirm this result and provide a more 
accurate estimate of the density contrast between the two 
components. 

In order to provide a further check on our results, we 
have measured the total [SII]AA6716,6731 fluxes by fitting 
only a single component to each line in the blend. This leads 
to ^|2if =1.24±0.06 and a density of rie = (1.9to;6) x 10* 
m~'^ (assuming re = 10,000) — consistent with Clark et al. 
(1997), who obtained n^ = (1.7±0.5) x 10* m"^. Note that, 
if the narrow component has a low density — as tentatively 
suggested by our multi-Gaussian fits to the blends — and 
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Kinematic components resolved in PKS2250-41 



* Spatially extended along the radio axis 
I and perpendicular to it. Not detected 

/ beyond the edges of the radio lobe. 

Broad ^ t-Flux and FWHM peak at the radio hot spot position 
Y ^ Large FWHM (500-900 km/s) 

* Low ionization level 

^Detected across the whole of the arc 
and far beyond the edge of the radio lobe 



Red // 
Narrow 



* Flux curve very similar along the radio axis, 
the 2"N and the 2.5 "S slit positions 

* FWHM~60-200 km/s, rather constant in space 
^Relatively high ionization level 

* Detected along the radio axis and in NW clump 
up to the edge of the radio lobe 



Blue 



^ * Relatively high ionization level 

Figure 6. Summary of the properties of tlie tliree kinematic components detected in the Western Arc of PKS2250-41. 



taking into account that the temperature of the broad com- 
ponent is higher, then Ue ~ 1.9 x 10* ni~^ represents a lower 
limit on the density of the broad component. 



4 DISCUSSION 

4.1 The nature of the different kinematic 
components 

The cartoons presented in Fig. 10 show 2-D pictures of the 
Western arc, with the spatial variation of the flux, FWHM 
and ionization level of the gas indicated. Darker squares in- 
dicate higher values. The relative intensity scale is approxi- 
mately correct for each cartoon, but not from cartoon to 
cartoon, (i.e. the same colour intensity in a given cartoon 
implies same value, but not in two different cartoons). The 
cartoons summarize the properties of the spatial variations 
of the different components in a more visual way. Notice 
the very similar line flux curves, ionization and kinematic 
properties of the red narrow component across the arc; the 
constancy of the flux and FWHM of the broad components 
across the arc, except at the position of the hot spot; and 
the presence of the blue component along the radio axis (and 
detached blob) with rather constant flux and FWHM. No- 
tice also the increase in the ionization level for all the com- 



ponents as the distance to the nucleus of the host galaxy 
decreases. 

The broad component 

The similarities in the properties of the broad compo- 
nent measured in all spatial positions, except in the proxi- 
mity of the hot spot, suggest that the broad components 
detected in all spectra have a common origin. This origin 
is clearly associated with the interaction between the radio- 
emitting structures and the ambient gas because: 

- the large velocity width is incompatible with pure 
gravitational motions at the position of the hot spot, and 
is also rather large at the other positions; 

- the very low ionization level is as expected from the 
compression effect of the shocks and/or the influence of 
shocks in the emission line processes; 

- the flux and line width peak at the position of the 
radio hot spot. Again, this is what we expect, since the hot 
spot is the working surface of the radio jet and the effects 
of the interaction are expected to be stronger here. 

- the high electron temperatures revealed by the [OIII] 
lines are consistent with the idea that the broad component 
is cooling gas behind a shock front. Although the density 
results are not conclusive, they suggest that the shocked gas 
has been compressed. 

- the faintness of the Hell line also suggests the effects 
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2250 PA270 6250 
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Wavelength Wavelength 

Figure 7. The fit to the H'y+[OIII]\4363 lines. Top panels: fit to the spectrum extracted from the radio axis sUt. Middle panels: fit to 
the spectrum extracted from the cross cut slit position. Bottom panels: fit to the spectrum sum of the two previous spectra. Left: Solid 
smooth line represents the fit to the data. Right: Solid smooth lines represent the individual components (narrow and broad). 
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Figure 9. The fit to the [SII]XX6716,6731 doublet. Left: Solid smooth line represents the fit to the data. Right: Solid smooth lines 
represent the individual components (narrow and broad). 
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Figure 8. The temperature: [OIII]A5007/[OIII]A4363 vs. pre- 
dicted T/10'' temperature. The results of the fits to the three 
spectra (see text) are presented. The broad component is repre- 
sented as an open circle and the narrow component as a solid 
triangle. The star represents the measurements from the total 
flux of the lines. All fits are consistent with the broad component 
being hotter. 



of shocks, since this property is predicted by shock models 
and similarly weak Hell lines are observed in other jet-cloud 
interaction targets {e.g. 3C171, Clark et al. 1998) 

The possibility that the broad component is scattered 
light is rejected since it is detected in the forbidden lines 



and the polarization level of the continuum is low (Dickson 
et al. 1995). 

The spatial extent of the broad component (~4.86 arc- 
seconds or ~30 kpc in the N-S direction) indicates that the 
effects of the interactions reach regions as far as ~15 kpc 
in the direction perpendicular to the radio axis — approxi- 
mately the width of the radio jet cocoon. 

The blue narrow component 

The apparent association of the blue narrow component 
with the radio axis along PAO, the abrupt disappearance of 
this component at the edge of the radio lobe, and the drop of 
the flux at the position of the hot spot suggest that this com- 
ponent may be affected by the interaction. There is marginal 
evidence from the [OIII](5007)/H/3 ratio along PAO (Fig. 3) 
that this component has a lower ionization state than the 
red narrow component. If confirmed, this result is consistent 
with the idea that the ionization of the blue narrow compo- 
nent has been affected by the compression and/or ionization 
effects of the jet-induced shocks. 

Alternatively, this component may simply represent a 
separate velocity system in the ambient gas of the host 
galaxy, with a higher density or different distance from the 
nucleus resulting in a lower ionization state. In this case, the 
apparent association with the radio axis along the cross-cut 
slit position would be coincidental. 

The red narrow component 

The emission of the red narrow component is characteri- 
zed by little variation in properties across the arc, including 
the radio axis and the position of the radio hot spot. The line 
widths are consistent with gravitational motions (~150 km 
s~^). The ionization level is high relative to the broad com- 
ponent. None of these properties present changes related to 
the radio structures. Moreover, this component is detected 
far beyond the radio hot spot (and beyond the edge of the 
radio lobe). 
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Figure 10. Cartoon showing tlie spatial variation of the flux (top panels), FWHM (middle horizontal panels) and ionization level (lower 
panels) for the different kinematic components detected in PKS2250-41: Red narrow component (left panels), broad component (middle 
vertical panels), blue narrow component (right panels). The position of the hot spot has been indicated with an 'x'. 



All of these properties suggest that the red narrow com- 
ponent is emitted by ambient material which has not yet 
been shocked. 

If the material emitting the red narrow component had 
passed through the shocks and subsequently cooled behind 
the shock fronts, we would expect a sudden decrease in the 
ionization state of this component at the position of the 
edge of the radio lobe, due to the compression effect of the 
shocks. As the H/3 line is too faint beyond the radio lobe, we 
have used the [OII]/[OIII] ratio to compare the ionization 
level beyond and behind the edge of the radio lobe. We have 
extracted a spectrum from a 1.6x1.5 arcsecond aperture (2 
pixels in the spatial direction) immediately behind the edge 
of the radio lobe and another spectrum from a 2.4x1.5 arc- 
second aperture (3 pixels) immediately beyond it. We obtain 
[OII]/[OIII]=2.1±0.2 behind and 2.0±0.3 beyond the edge 
of the radio lobe. Therefore, there is not the sudden change 



in the ionization level that we would expect if the red narrow 
component represents cooled post-shock gas. 

Interestingly, the spatial distribution of this component 
— particularly the evidence for edge brightening along the 
slit position perpendicular to the radio axis (Fig. 5) — shows 
that it is associated with the arc structure observed in the 
[OIII] emission line image (see Fig. 1). The question then 
arises as to why this component has an arc shape — remi- 
niscent of a bow shock — yet there is no spectroscopic ev- 
idence that the gas emitting the red narrow lines has been 
disturbed by the shocks. We will return to this question in 
section 4.4.1. 

Therefore, we have resolved kinematically the emission 
from the shocked gas and the ambient gas in PKS2250-41- 
Both gaseous components present different kinematic and 
physical properties. 
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We note striking similarities between PKS2250-41 and 
two other jet-cloud interactions candidates which have re- 
cently been studied in depth: 3C171 (Clark et al. 1998) 
and PKS1932-46 (Villar-Martin et al. 1998). In the case of 
3C171, broad and narrow components are observed across 
the extent of the interaction and, just as we observe in 
PKS2250-41, the broad component has a low ionization state 
while the narrow component has a high ionization state and 
extends beyond the radio structure. The HeIIA4686 line (in- 
tegrated flux) is also weak. Although the case for a strong 
jet-cloud interaction in PKS1932-46 is weaker, we also find 
evidence in that object for a narrow, high ionization com- 
ponent which extends well beyond the radio lobe, and a 
broader, lower ionization component which may be related 
to a jet-cloud interaction. Remarkably, this similarity in 
the kinematic/ionization properties persists despite a wide 
range in the emission line morphologies of the individual 
sources. 



AGN ILLUMINATION 



Photoionized ambient gas 



^ \ ^ 




Bow shock 



4.2 The scenario 

We now propose a scenario to explain the origin and pro- 
perties of the different kinematic components resolved in 
PKS2250-41. This scenario could be applied to other radio 
galaxies in which jet-cloud interactions are taking place. 

Hydrodynamical models of an advancing jet in a 
gaseous medium show that a bow shock is created that ex- 
pands transversely to the jet axis {e.g. Hartigan et al. 1987, 
Taylor et al. 1992) (see Fig. 11 based on Bicknell et al. 1997). 

The models also predict that the shocked gas is first 
heated to temperatures beyond lO^K and cools slowly in 
a phase ruled by the decreasing pressure (e.g. Ferruit et 
al. 1997). When the gas cools to a temperature of a few 
lO^K a catastrophic cooling occurs. The temperature drops 
very quickly, while the density increases. When the gas cools 
down to a few x lO^'K we will be able to observe the line emis- 
sion. The post-shock gas contains a significant range of gas 
temperatures and the [OIIIJA4363 line luminosity is domi- 
nated by emission from the 20,000-30,000 K zone. This is 
consistent with the high observed temperatures in the broad 
component and the possible compression suggested by the 
low ionization state and the density measurements. 

The gas passing through the shock front will be acce- 
lerated and its kinematics will be perturbed. The effects are 
likely to be stronger at the position of the hot spot, the 
working surface of the radio jet. This is also what we observe 
(large FWHM with maximum value at the position of the 
hot spot). 

Therefore, the properties of the broad component are 
consistent with it representing gas which has passed through 
the shock front and been accelerated, heated and com- 
pressed. Our results suggest that we are observing the gas 
while the dramatic cooling is happening (it is hotter than 
the ambient gas). 

The blue narrow component is reminiscent of linear 
structures associated with the jets in some Seyfert galaxies. 
Such structures are characterized by high ionization level 
compared to the gaseous cocoon around the radio lobes 
(Capetti et al, 1996). The similarities with the properties of 
the blue narrow component in PKS2250-41 suggest a simi- 
lar origin, although the nature of such linear features is not 
understood with any certainty. 



Bow shock advance 



Figure 11. The bow shock model. (Bicknell et al.) 



On the other hand, the properties of the red narrow 
component are consistent with the idea that it is emitted 
by warm gas which has not yet been shocked. It reveals the 
properties of the underlying quiescent gas. However, uncer- 
tainties remain concerning how this component is ionized 
(see §4.4.1). 



4.3 The acceleration mechanism 

We have measured line widths of up to ~900 km s^^ 
(FWHM) in PKS2250-41. The real velocities could be much 
larger, since, given the relatively weak radio core (Clark et 
al. 1997), it is likely that the jet is moving close to the plane 
of the sky, and any systematic motions will be induced in 
that plane. Line widths of up to several thousands km s~^ 
have also been measured on the EELR of many high redshift 
radio galaxies {e.g. McCarthy et al. 1996, Rottgering et al. 
1997). 

One of the big challenges for models of interactions be- 
tween the radio and optical structures is to explain how the 
warm gas is accelerated to such large velocities. The jet prop- 
agates through a two phase medium, with warm (T'^IO^ K) 
clouds embedded in a hot (T~10^ K) ambient phase, thus 
corresponding to a density contrast of ~ 1000 if the clouds 
are in pressure balance with the hot phase. There are several 
possibilities for accelerating the warm gas. These include the 
following: 

• Interaction between the cloud and the bowshock 
ahead of the jet working surface: in the absence of en- 
trainment processes (see below) the cores of the warm clouds 
will be accelerated to Vw ~ Vs{ph/pw)'^'^ (e.g. Klein et al. 
1994) where, Vw is the final velocity of the warm clouds, Vs 
is the shock speed in the hot phase, ph is the preshock den- 
sity of the hot phase, and p^ is the preshock density of the 
warm phase. Estimates of hot spot advance speed (0.01c- 
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Figure 12. Velocity distributions of material along the line of sight for a cloud (lOpc diameter) that is impacted by an expanding 
bowshock associated with the radio plasma. Left and right panels correspond to two stages during the disruption (t = 2 X 10* yr on the 
left and 4 X 10'* yr on the right). The histograms show the total amount of mass, or relative column densities, of all the material with 
temperature T < lO'' K along the line of sight in each velocity bin, with three inclinations of the shock velocity vector to the line of 
sight (i represents the orientation of the entrained cloud material with respect to the plane of the sky). The parameters of the cloud, the 
ambient medium and the bowshock are appropriate for a rich group of galaxies (consistent with the X-ray properties of PKS2250-41: 
Clark et al. 1997); an initially homogeneous, spherical cloud of temperature T ~ 15, OOO-R', in pressure balance with an ambient medium 
of T ~ lO'^K and pressure P ~ 4 X 10~^^ dyn cm~^, is impacted by a shock moving at 5000 km s~^. The simulation was carried out 
using the Zeus-3D code, in a 3-dimensional volume of 192x64x 64 cells, with the cloud diameter initially resolved by 24 cells. 



O.lc, e.g. Scheuer 1995) lead to shock velocities in the range 
3x10'^ < Vs <3xl0* km ~^ . For warm clouds in pressure 
balance with a hot halo (T~10^ K) the density contrast is 
~1000 between the hot and the warm phases. Therefore, the 
range of velocities of the accelerated warm clouds is 100-1000 
km ~^ . However, noting that this is the maximum we would 
expect for the working surface of the jet, that much of this 
velocity is likely to be directed in the plane of the sky (the 
direction of the jet), and that the cocoon away from the 
hotspot is likely to expand much more slowly than at the 
working surface, it is unlikely that the bowshock accelera- 
tion is sufHcient to explain the large line widths observed in 
PKS2250-41 and other jet-cloud interaction candidates. 

• Entrainment in the hot post-shock wind. An in- 
teresting possibility is that the clouds are entrained in the 
hot wind behind the shock. Instabilities will arise at the 
interface between the hot, fast flowing gas and the warm, 
much slower moving clouds. These instabilities will lead to 
entrainment and acceleration of warm material in a hot wind 
(e.g. Stone & Norman 1992; Klein et al. 1994; Dai & Wood- 
ward 1994). Koekemoer (1998) has recently completed a grid 
of models in the specific context of clouds with high density 
contrast (~ 1000) in the hot gas around radio galaxies, which 
are impacted by shocks with properties typical of expanding 
radio cocoons. In Fig. 12 we show results from one such sim- 



ulation carried out in the specific context of PKS 2250-41: 
a cloud initially at T ~ 15, QQQK embedded in a hot inter- 
stellar medium with conditions appropriate to a rich group 
of galaxies is impacted by a shock moving at 5000 km s^^, 
expected to be typical of the bowshock expansion velocity. 
The velocity profiles of the stripped gas (with T < 10® K) 
are shown. Although this is a relatively simple scenario, not 
yet taking into account time averaged profiles or different 
distributions of cloud sized along the line of sight, it is nev- 
ertheless clear than the velocities of order ~ 1000 km s^^ 
can be produced in this way. By integrating across the sur- 
face of the shock, for a range of cloud properties and stages 
of interaction, it would clearly be feasible to produce broad 
lines similar to those observed in PKS2250-41. 

• Entrainment of clouds in the turbulent boun- 
dary layers between the radio jet and the ISM or 
the jet cocoon and the ISM (e.g. Sutherland et al. 1993). 
This mechanism could also produce large line widths and 
accelerate the clouds to a significant fraction of the jet speed. 

Therefore it is likely that, in addition to the initial shock 
acceleration, an entrainment process is required to explain 
the large widths of the broad emission lines. Higher spatial 
resolution optical and radio observations will be required 
to determine whether the gas is being entrained in the hot 
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Figure 13. Diagnostic diagrams: AGN photoionization models -> Solid line {logU values are indicated in the logfNIIJ/Ho vs. 
lg([OI]/[OIII]) diagram). Shock models (cooling gas) -> Long dashed lines. Shock+precursor models — > dotted lines (values of the shock 
velocity Vg and the magnetic parameter B are indicated in the log(HeII/H/3 vs. lg([OI]/[OIII]) diagram). Nucleus 

— > Solid square. Broad component -> open circle. Narrow component -> open triangle. 



post-shock wind or in the turbulent boundary layers of the 
radio jets. 



4.4 Ionization mechanism: 
photoionization 



shocks vs. AGN 



Shocks can influence the ionization of the gas in two ways 
(Sutherland ei al. 1993, Dopita & Sutherland 1996): 1) the 
initial heating of the gas as it passes through the shock front, 
leading to coUisional ionization of the gas and, eventually, 
to forbidden line emission as the gas cools radiatively; 2) 
the generation of a strong UV continuum in the hot post- 
shock gas which can ionize both the precursor gas and also 
the gas which has cooled radiatively behind the shock front. 
However, the galaxy also contains an active nucleus which 
is likely to generate a hard continuum capable of ionizing 



the gas. Which mechanism dominates the ionization of the 
extended emission line gas in PKS2250-41: shocks or AGN 
photoionization? 

Taken at face value, the high electron temperature mea- 
sured for the broad component in the western arc provides 
strong evidence that the gas emitting the broad component 
is cooling behind a shock front. Such high temperatures are 
predicted by models for shocked, cooling gas, but are difficult 
to explain in terms of conventional, radiation bounded AGN 
photoionization models. Although photoionization models 
which include matter-bounded clouds, or clouds photoioni- 
zed by the hard continuum which has been partially ab- 
sorbed by a matter-bounded component, can produce high 
electron temperatures (Binette et al. 1997), such models can- 
not explain a temperature as high as Te > 20, OOOK and a 
low HeII(4686)/H/3 ratio simultaneously. 
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In order to check the consistency of our results with the 
ionization models in greater depth, we have compared the 
measured values for the line ratios of the different gaseous 
components in PKS2250-41 with the prediction of both AGN 
photoionization and shock models. The results are shown 
in the diagnostic diagrams presented in Fig. 13. We have 
used the multi-purpose photoionization code MAPPINGS I 
(Binette et al. 1993a,b) to generate the AGN model predic- 
tions. We have assumed a density of 100 cm~'^ at the illumi- 
nated face of the clouds, solar abundances, and a power-law 
spectral index for the photoionizing continuum of a = 1.5 
{Fi, oc 1^'°")- The results are presented as continuous lines in 
Fig. 13 as sequences in the ionization parameter U, defined 
by 



U 



1 

cnn 



hv 



-du 



where c is the speed of light, uh is the density of the 
gas in the front layer, vq is the Lyman limit frequency, and 
(fiu is the monochromatic ionizing energy flux impinging on 
the slab. Note that this same sequence is able to reproduce 
the general trend defined by the optical line ratios of low 
redshift radio galaxies (Robinson et al. 1987). 

For the shock models we have used the published results 
of Dopita and Sutherland (1995), whose models include both 
the gas cooling behind the shock front, and the gas pho- 
toionized by the hot post-shock gas (the photoionized pre- 
cursor). The two main parameters which influence the pre- 
dicted spectrum are the velocity of the shock [vg C [150,500] 
km s~^) and the magnetic parameter (B C[0,4] ^G cm~"^"). 
Adopted density is n(H) = lcm~"^ and solar abundances. The 
spectrum predicted for the cooling gas is presented in the 
diagnostic diagrams in Fig. 13 as long dashed lines. Each se- 
quence corresponds to a flxed B value and changing Vs . The 
spectrum predicted for the precursor is presented as dotted 
lines {vs C[150,500] km s"^ and fixed B =1 /xG cm'^^^). 

Fig. 13 shows that the broad and narrow component oc- 
cupy very different positions in the diagnostic diagrams, due 
mainly to a lower ionization level of the broad component. 
Although the line ratios for the broad component are repro- 
duced by the AGN photoionization models on many of the 
diagrams, an important problem of AGN models is the failu- 
re to explain the high electron temperatures and the faint- 
ness of the HeIIA(4686) line for the broad component. We 
plot in Fig. 14 a new diagnostic diagram: [OIII]A5007/4363 
vs. IIeII/II/3. This diagram clearly distinguishes the narrow 
and the broad components. The AGN photoionization mod- 
els do not appear in this diagram at all because they predict 
high [OIII] 5007/4363 ratios {i.e. corresponding to low elec- 
tron temperatures) , which fall above the top of the diagram. 
Taking Figures 13 and 14 together it is clear that the line 
ratios of the broad component are entirely consistent with 
the idea that this component represents collisionally ioni- 
zed gas which is cooling behind a shock front; the shock 
velocities able to reproduce the position of the broad com- 
ponent in the diagnostic diagrams are in the range 200-300 
km s~^. The AGN photoionization models fail to produce 
such a good overall fit to the emission line spectrum. Thus 
we have strong evidence that the jet-induced shocks not only 
compress, but also ionize, the warm clouds (i.e. the jets have 
a significant energetic input). 
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Figure 14. The temperature: [OIII] A5007/ [OIII] A4363 vs. 
HcII/H/3. The broad (open circle) and the narrow (solid trian- 
gle) components are clearly separated in this diagram. Dotted 
lines: shock -fprecursor models. Dashed lines: shock models. AGN 
photoionization models lie outside the diagram. They produce too 
high [OIII] 5007/4363. 



In contrast to the situation for the broad components, 
the AGN and shock precursor photoionization models pro- 
duce good fits to the line ratios of the narrow components 
on most of the diagrams. Any remaining discrepancies be- 
tween the measured and predicted [OIII] 5007/4363 ratios for 
the narrow components can be resolved by adding a matter- 
bounded component (Binette et al. 1997). 



4.4-1 The ionization of the red narrow component 

We noted in section 4.1 that there is an apparent inconsis- 
tency between, on the one hand the association of the red 
narrow component with the arc structure, and on the other 
the lack of clear spectroscopic evidence that the material 
emitting this component has been shocked. One possible ex- 
planation — consistent with the emission line spectrum of 
this component — is that the red narrow component repre- 
sents shock-photoionized precursor gas. 

The shock photoionization process will be most efficient 
when the precursor gas is situated in the clouds that are be- 
ing shocked, since the precursor gas will then absorb a large 
fraction of the ionizing photons generated in the shocks. As 
the fast shocks being driven through the hot ISM by the ra- 
dio source engulf denser clumps of warm gas, slower shocks 
will be driven through the warm clouds (see section 4.3). 
The X-ray photons generated in the hot gas immediately 
behind the shocks in the warm clouds will not be efficient 
at photionizing the warm gas, however, provided that the 
gas cooling time {tcooi) is short compared with the shock 
crossing time {tcroas), a cooling zone will form downstream 
of the shock which will emit lower energy photons, capa- 
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ble of ionizing the precursor gas efficiently. Thus we require 
tcross > tcooi for efficient precursor ionization. Using the 
approximate formula for the cooling time given in section 
2.3 of Klein et al. (1994) this leads to the following cons- 
traint on the diameter [D) of warm clouds: 



->-<wn 



v„ 



200 km s- 



pc 



where no is the precursor density (in units of nC^), and Vw is 
the velocity of the slow shock through the warm cloud. Thus, 
for a precursor density typical of warm clouds in pressure 
equilibrium with the hot X-ray halo of a group of galaxies 
(no ~ 1 X lO^rn"^) and a shock speed consistent with the 
measured emission line ratios of the broad component (Vw ^ 
200 km s~^), we require D > 1.8 pc for efficient precursor 
photoionization. Although we have no direct indication of 
the diameters of the warm clouds interacting with the radio 
components in PKS2250-41, it is notable that typical spiral 
galaxies contain populations of molecular cloud complexes 
with diameters that are orders of magnitude larger than this 
limit [e.g. Mihalas & Binney 1981, Chapter 9). Therefore, 
it is likely that the PKS2250-41 system will contain at least 
some clouds for which the precursor photoionization process 
envisaged above is feasible. 

We now estimate the mass flow rate required to main- 
tain the narrow HP luminosity (Lhib)- For a strong shock, 
the H/3 luminosity of the photoionized precursor gas is re- 
lated to the shock velocity through the warm clouds (Vw) 
and the mass flow rate through the shock (M) by the fol- 
lowing equation (adapted from equation 4.4 of Dopita & 
Sutherland 1996): 



Lhi3 = 1.45 X 10^^ [ 
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For our chosen cosmology, we estimate a total H/3 lu- 
minosity for the red narrow component emitted by the arc 
of Lhp = 2.2x10^"* W . Substituting this into the expression 
for the HP luminosity we obtain: 

^ = ^-^><^°l200l^)"' ^«^^"- 

To put this mass flow rate into context, if the radio 
source requires ~10® yr to traverse the region of the arc, and 
Vs ^ 200 kms~^, the total amount of material "consumed" 
in the shock will be 1.5 x 10^ M© — comparable with the 
total HI content of a typical spiral galaxy. Although large, 
this is consistent with the idea that the phenomena we ob- 
serve in PKS2250-41 are the result of a direct interaction 
between the radio-emitting components and a companion 
galaxy. Note also that there is more than enough mechani- 
cal energy in the jets to power the emission line regions (see 
Clark et al. 1997). 

As well as the red narrow component directly associated 
with the arc, we also see a more extended red narrow com- 
ponent which extends well beyond the western radio lobe. 
Given its spatial extent, this component is unlikely to be 
emitted by clouds which have been engulfed by the fast 
shock driven through the hot ISM by the radio components. 
Nonetheless, it remains possible that this component is pho- 
toionized by photons from the warm, shocked clouds behind 
the shock front, provided that the covering factor of the 
extended component is large enough to intercept a signifl- 



cant fraction of ionizing photons generated by the shocks. 
However, a potential problem with this scenario is that the 
effective ionization parameter of the shock continuum will 
be much larger for the warm clouds in the arc than for the 
more extended warm clouds, yet the measured ionization 
state for these two components is similar (see section 4.1). 
While it is possible to get around this problem by assuming 
that the more extended narrow component is associated a 
lower density ISM than the narrow component emitted by 
the arc, it seems unlikely that the density gradient in the 
ambient ISM would exactly compensate for the change in 
the flux density of the shock continuum with distance from 
the shock front. 

In view of this potential problem, we should also con- 
sider the alternative possibility: that the emission line arc 
represents the intrinsic structure of the AGN-photoionized 
ambient gas i.e. the observed arc-like morphology is not a 
consequence of the shock, and the resemblance to a bow 
shock is misleading. Arc-like structures in the ISM are read- 
ily produced as a consequence of galaxy mergers and in- 
teractions, and PKS2250-41 shows morphological evidence 
for interactions with at least one companion galaxy (Clark 
et al. 1997). However, it would seem a remarkable coinci- 
dence if, purely by chance, such a merger remnant exactly 
circumscribed a radio lobe. Furthermore, arc-like structures 
circumscribing the radio lobes have been observed in other 
radio galaxies with jet-cloud interactions (e.g. PKS 1932-464, 
Villar-Martm et al. 1998 and Coma A, Tadhunter et al. in 
preparation). 

Thus, no explanation for the ionization and distribution 
of the red narrow component seems entirely satisfactory, and 
this issue clearly warrants further investigation. 

4.5 Shocked and ambient gas in high redshift 
radio galaxies 

PKS2250-41 shares many similarities with high redshift ra- 
dio galaxies. Therefore, what we learn studying this ob- 
ject will help us to understand the observed properties of 
very distant radio galaxies. Like many high-z radio galaxies, 
PKS2250-41 shows perturbed kinematics and a close correla- 
tion between the EELR and the radio structures. McCarthy 
& van Breugel (1989) noted that the brightest EELR in radio 
galaxies at z>0.7 tend to be associated with the radio lobe 
closest to the nucleus of the host galaxy which moreover, 
is systematically more depolarized (Pedelty et al. 1989, Liu 
& Pooley 1991). PKS2250-41 also shows these asymmetries 
(Clark et al. 1997) which may arise, as proposed for high-z 
radio galaxies, as a consequence of an asymmetric gaseous 
environment surrounding the host galaxy. 

We have shown that the strong interactions between 
the radio jet and the ambient gas imprint special features in 
the observed properties of the ionized gas, which differ from 
the ambient (non shocked) gas properties. Do we see these 
features in HzRG?. 

At least some HzRG show, in addition to the high ve- 
locity gas associated with the radio structures, large halos 
of ionized gas extending far beyond the radio structures. 
Interestingly, the kinematics of this gas are rather quies- 
cent compared to the inner structures. Examples include 
1243+036 (2 =3.6, van Ojik et al. 1996), 3C368 [z =1.1, 
Stockton et al. 1996) and 4C41.17 {z =3.8, Dey et al. 1997, 
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Chambers et al. 1990). The high velocity gas in both 3C368 
and 4C41.17 presents low ionization level compared to the 
low velocity gas. 

1243+036, 3C368, 4C41.17 are examples of distant radio 
galaxies which show the presence of two types of ionized 
gas. The similarities with PKS2250-41 suggest that the na- 
rrow line gas (at least the gas detected beyond the radio 
structures) is probably quiescent ambient gas which has not 
interacted with the radio structures and is ionized by the 
AGN continuum and/or the shock UV continuum. On the 
other hand, the high velocity gas is shocked gas which has 
been heated, compressed and kinematically disturbed by the 
shocks. 

Therefore it seems likely that the ambient and shocked 
gas have also been resolved in some of the high-z objects, and 
they present differences in their kinematics and ionization 
properties consistent with the results obtained for PKS2250- 
41. A similar kinematic study as described in this paper will 
be very fruitful in these objects. It will allow us to isolate 
the non-shocked ambient gas and study the properties of the 
gaseous halos in elliptical galaxies in the young Universe 
as well as to understand the way the observed properties 
are determined by the interaction with the radio-emitting 
plasma. 



the shock front. The arc shape circumscribing the radio lobe 
suggests that the second mechanism is dominant, although 
a large mass flow rate through the shock is required. 

The similarities between the kinematic properties of 
PKS2250-41 and some high-z radio galaxies suggest that 
the ambient and the shocked gas have also been resolved 
in the more distant sources, while the similarities between 
PKS2250-41 and other jet-cloud interaction radio galaxies 
(3C171, PKS1932-464) in spite of their very different mor- 
phology, suggest that jet-cloud interactions have rather con- 
stant effects on the physical and kinematic properties of the 
emission-line gas in powerful radio galaxies. 
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5 SUMMARY AND CONCLUSIONS 

We have studied the effects of the interaction between the 
ambient gas and the radio-emitting plasma in PKS2250-41. 
The results of this research provide clues for understanding 
the processes which determine the observed properties of the 
majority of very distant radio galaxies and to determine the 
role played by the interactions between the radio and optical 
structures in active galaxies in general. 

We have resolved kinematically the emission from the 
gas which has passed through the shocks induced by the 
radio plasma and the emission from the ambient precursor 
gas. The properties of the two main gaseous components 
show marked differences. 

Our analysis shows that the effects of the interactions 
can spread at least 15 kpc perpendicular to the radio axis 
— across the full extent of the radio lobes. The interac- 
ting gas presents broad line widths (FWHM>900 km s~^), 
a low ionization level and weak HeIIA4686 compared to the 
ambient gas. We have presented the first evidence that the 
gas emitting the broad component is hotter (T ~30,000 K) 
than the ambient gas (T ^^15,000 K). All of these properties, 
together with the morphological association with the radio 
structure, are consistent with the broad component being 
emitted by gas which has been heated and coUisionally ioni- 
zed by a jet-induced shock and is now radiatively cooling 
behind the shock front. This represents compelling evidence 
that jet-cloud interactions not only compress and accelerate 
the interacting warm clouds, but also ionize them. 

The ambient gas is characterized by a relatively high 
ionization level and a more quiescent kinematics, consistent 
with gravitational motions. The properties of this ambient 
gas can be explained in terms of pure gravitational motions 
and photoionization by a hard continuum source. It is not 
yet clear whether this component is photoionized by the 
AGN or by UV photons emitted by the gas cooling behind 
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